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Abstract 

A new concept of remote detection of concealed radioactive materials has been recently 
proposed PQ-[3]. It is based on the breakdown in air at the focal point of a high-power beam 
of electromagnetic waves produced by a THz gyrotron. To initiate the avalanche breakdown, 
seed free electrons should be present in this focal region during the electromagnetic pulse. This 
paper is devoted to the analysis of production of free electrons by gamma rays leaking from 
radioactive materials. Within a hundred meters from the radiation source, the fluctuating free 
electrons appear with the rate that may exceed significantly the natural background ionization 
rate. During the gyrotron pulse of about 10 microsecond length, such electrons may seed the 
electric breakdown and create sufficiently dense plasma at the focal region to be detected as an 
unambiguous effect of the concealed radioactive material. 

1 Introduction 

A new concept of remote detection of concealed radioactive materials has been recently proposed 
PQ. This concept is based on breakdown in free air at the focal point of a high-power beam of 
electromagnetic waves. To realize such breakdown, the wave amplitude in a focal region should 
exceed the breakdown threshold. To initiate the avalanche breakdown process, seed free electrons 
should be present in this focal region during the electromagnetic pulse. When the wavelength is 
short enough (sub-THz or THz frequency range) the wave beam can be focused in a spot with 
dimensions on the order of a wavelength and, then, the total volume where the wave amplitude 
exceeds the breakdown threshold can be rather small. This fact allows one to realize conditions pQ 
when the breakdown rate in the case of the ambient electron density is rather low. So in the cases 
when the breakdown rate is much higher than the expected one can conclude that in the vicinity of 
a focused wave beam there is a hidden source of radioactive material which ionizes the air. Some 
issues important for realizing this concept are discussed in Refs. [2] and [3]. 

In the present paper, we analyze another issue important for realizing this concept which was 
only briefly discussed in previous references, viz. propagation of 7-rays in air and production of 
free electrons by these quanta, mainly due to inelastic Compton scattering. It should be noted 
that, in principle, due to the atomic reaction radioactive materials emit both MeV-scale 7-rays 
and /3-electrons. For example, a disintegrating atom of 27C0 produces one 2.505 MeV /3-electron 
and two 7-quanta with the energies 1.173 and 1.332 MeV [5]. Metal shielding and container walls 
largely stop /3-electrons, while absorbing only a fraction of 7-rays. All unabsorbed 7-rays, even 
scattered ones, will eventually leak through the walls and propagate in air. An example of energy 
distribution of 7-rays and electrons outside the container wall is shown in Fig. [T] reproduced from 
Ref. [2]. 
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Figure 1: Results of Monte-Carlo simulations (MCNP code) for 1 Curie of 27C0 showing the energy 
spectrum of electrons (red) and photons (blue) in a 1cm 3 tally volume outside a 0.5 cm thick steel 
shipping container. Solid lines are the results of calculations for the case when electron cutoff energy 
inside the steel wall adjacent to tally volume is set to 1 keV, while dotted lines are for electron 
cutoff energy set to 100 keV. 



The energetic 7-quanta escaped from the container will occasionally collide with air molecules, 
experiencing Compton scattering and photoelectric absorption, as discussed below. Their propa- 
gation distance can be estimated based on the interaction cross-section a given, e.g., in Ref. [3J on 
Figures 2.50 (d) and (e) for nitrogen and oxygen, respectively. (One of these very similar figures 
was reproduced in Ref. [2].) From these figures one can draw, at least, two conclusions. First, for 
energies between 0.1 MeV and 1.2 MeV shown in Fig. 1, the dominant effect in the production of 
free electrons is the Compton scattering (this conclusion agrees with Figures 2.42, 2.43 of Ref. [1]). 
Second, the interaction cross-section a, as the function of energy, gradually increases from about 
0.055 cm 2 /g at 1.2 MeV to about 0.16 cm 2 /g at 0.1 MeV. Correspondingly, the mean free pass 
of 7-rays in air (n a = 2n m ~ 6 • 10 19 atoms/cm 3 ) defined as I = l/(n a a) varies from ~120 m for 
1.2 MeV quanta down to about 40 m for 100 keV quanta. When the energy of 7-quanta decreases 
down to 30 keV and below, the photoelectric absorption becomes the dominant effect. 

It is important that both the Compton scattering and photoelectric absorption ionize air 
molecules, producing high-energy electrons within the MeV range. These primary electrons, which 
should not be confused with the /3-electrons directly from radioactivity, cannot propagate too far: 
they will be absorbed by air within a few meters [3J. During their short lifetime, however, most 
of these primary electrons produce thousands of secondary ones. It is this process that determines 
the local rate of total electron production required for the avalanche breakdown [3]. As discussed 
above, all unabsorbed 7-quanta propagate away from the radioactive source, continuing to produce 
free electrons. Hence to estimate the total rate of free-electron production at various distances from 
the radioactive source, one needs to analyze the collisional propagation of 7-quanta in air. In this 
paper, we analyze the photon propagation in the near-source zone, i.e., at distances less or of order 
of the mean free pass I estimated above. 
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Figure 2: Kinetic energy, KN differential cross section and angular differential cross-section of 
gamma photons. 
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2 Production of free electrons 



As propagating 7-quanta collide with air molecules, producing high-energy primary electrons, they 
experience a strong angular scatter and energy losses. The major electron production takes place 
when the photon energy £~ remains between 30 keV and 1.3 MeV, i.e., in the energy range where 
the Compton scattering dominates. In this range, one can successfully apply the Klein-Nishina 
(K-N) theory that describes Compton scattering of photons by free electrons [4]. Indeed, after a 
typical collision between the incident photon and an electron bound within an O2 or N2 molecule, 
the originally bound electron becomes free and acquires a significant momentum (energy). If the 
acquired energy £ 7 greatly exceeds the electron binding energy then both the scattered photons 
and released electrons will have momentum distributions fairly close to those if the originally bound 
electron was free and at rest. Furthermore, if the incident-photon energy is well above the maximum 
electron binding energy (the K-shell edge), Sk — 500 eV for the air molecules, then the K-N theory 
should be applicable to all bound electrons with the exception of small-angle scattering when the 
released electrons acquire a small energy, £ e < Ek [H Section 2.6.2.1]. 

As described in numerous textbooks, see, e.g., Ref. [U Section 2.6], in the frameworks of the 

K-N theory, the kinetic energy, T' = c [p 2 + m^c 2 ) — m e c 2 = £ 7 — £ 7 , of the outgoing electron 
and the scalar momentum of the scattered photon, p'^ = jp^j = £!y/c (for the incident photon with 
p 7 = |p 7 | = Sq/c) are given by 
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while the K-N differential scattering cross-section is given by 
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In (jTJ) — , 7 is the scattering angle, i.e., the angle between the momenta of the incident and 
scattered photons, p 7 and p^, which is rigidly related to p 7 and j? 7 , 

m e c m e c 

cos 7 = 1 ; — I , (3) 

fly = (0 7 , $ 7 ) is the solid scattering angle, m e is the electron rest mass, c is the vacuum speed of 
light, a = £ 7 /(m e c 2 ) = p 7 /(m e c) is the energy of the incident photon normalized to the electron 
rest energy and ro ~ 2.818 x 10~ 13 cm is the classical electron radius. 

The K-N formula ([2|) allows one to calculate the total cross-section and, therefore, the production 
rate of the primary electrons, 2Zp L \ r n m ca\<^, where n m = n a /2 ~ 3 • 10 19 cm -3 is the total density 
of air molecules and Zj^ir is the mean air atomic number (for the standard air composition of 79% 
N2 with = 7 and 21% O2 with Zq> = 8, we have Zp^ ~ 7.21). Of prime interest to us is the 
production rate, dN e /dt, of secondary (knock-on) electrons which will inevitably end up attaching 
to neutral molecules, i.e., forming negative ions [3]. Each primary electron with the kinetic energy 
T produces on average T/(34eV) secondary electrons [5]-[7]. Then, to determine the total average 
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rate of secondary electron production by photons of a given initial energy a, we should first integrate 
the K-N differential cross-section over all solid angles with the weight factor T' given by ([Taj) and 
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and then divide the result by 34 eV. The exact calculation of the last integral of yields a function 
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shown in Fig. [3] by the black solid line. Note that the function K(a) is proportional to the average 
"absorption cross section" [9j Eq. (8e-28)], 



-KrlK(a) 



a 



(5) 



In the energy range of interest, £ 7 < 1.332 MeV, i.e., a < a max w 2.61, the complex function K{a) 
can be quite accurately approximated by the dependence 0.37a (a a = 0.377rrQ w 9.23 x 10 -26 cm 2 ) 
shown in Fig. [3] by the red dashed line. This approximation allows us to describe the total average 
rate of secondary electron production by a simple formula: 
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where -F 7 (a) is the photon energy distribution normalized to the total photon density, n 7 = 
f F^(a)da, and expressed in terms of the full momentum distribution function / 7 (p 7 ) (Sect. [3]) 
as _F 7 (a) = m 3 c 3 a 2 J / 7 (p 7 , 7 )<ifi 7 . Thus, to evaluate the average rate of secondary electron pro- 
duction as a function of the distance from the radiation source, one needs to determine the photon 
energy distribution along with its distance variations. This requres a collisional kinetic theory for 
propagating photons. 



3 Photon kinetics in spherically symmetric approximation 

Radiated and scattered 7-quanta can be treated as a gas of incoherent and unpolarized photons. 
For photon energies above 30 keV, this gas is prone mainly to Compton scattering, while at lower 
energies the photoelectric absorption starts dominating. Both processes result in air molecule 
ionization, i.e., production of free electrons. The coherent Rayleigh scattering is at least an order 
of magnitude smaller and can be neglected. Also, the process inverse to photon absorption, viz., 
photon production during collisions of free electrons with air molecules, in the range of electron 
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energies we are dealing with (< 1.3 MeV), is negligibly small. Therefore, we will assume the spatial 
and momentum distribution of photons to be independent of any feedback from the free electrons. 

We will treat individual 7-quanta kinetically as classical particles localized in both real and 
momentum spaces. The photon momentum distribution, / 7 (p 7 ), may have a temporal (t) and 
spatial (r) dependences. Since photons propagate and scatter practically instantaneously we will 
ingore any temporal dependence. Being concerned mainly with distances from the radioactive 
material much longer than the container size a, we will treat this material as a point-like source of 
7-quanta. Ignoring possible effects of the ground, air inhomogeneity and container anisotropy, we 
will assume spherically symmetric propagation of photons. 

In the spherically symmetric geometry, the photon momentum distribution depends on one 
spatial variable, viz., the distance r of a localized photon from the point-like source, and two 
momentum components: the absolute value of the momentum, p 7 , and the momentum polar angle 
with respect to the radius vector at the given location, # 7 (the latter should not be confused with the 
angle between the scattered and incident photons, 7 ). Thus the photon momentum distribution 
is a function of three variables, / 7 (r,p 7 , # 7 ). Since the angle # 7 is not invariant with respect to the 
collisionless motion, it is more convenient to pass from 9~ f to the "impact parameter" p, 

( P 2 V /2 

p = rsin# 7 , cos # 7 (p) = ± ( 1 ^ J > (7) 

equal to the shortest distance from the radioactive source (r = 0) to the straight-line trajectory of 
a freely propagating photon. The parameter p is proportional to the photon angular momentum 
and remains invariant during the free photon propagation. It can be changed only during brief 
collisions that lead to Compton scattering. For any given positive p, we distinguish between two 
groups of photons with different signs of cosO-y(p) in ([7]): the "outgoing" photons with the '+' sign 
(corresponding to dr/dt > 0) and "ingoing" photons with the ' — ' sign (dr/dt < 0), as illustrated 
by Figure [H We will denote the two corresponding distribution functions by /if . 

In terms of /if (r,p 7 ,p), the total photon number density, n 7 , and radial flux density, g 7 (r), 
are given by 

2vr ^ f°° 2 [ r f±(r,p 1 ,p)pdp 
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The momentum distribution function of freely propagating 7-quanta with occasional spiky col- 
lisions satisfies the Boltzmann kinetic equation derived in Appendix lAl 

n 2 \ 1/2 df± 



±c(l-^J -£ r = L r f! r ± -v( Pj )f±, (9) 



where the left-hand side describes the collisionless propagation of photons, while the right-hand 
side (RHS) describes their collisions with air molecules. 

The first term in the RHS of ([9]) describes the kinetic arrival of photons to the elementary phase 
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volume d 3 p~f around given p 7 due to collisional scattering of incident photons from d 3 p' around p 
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where 



5 (cos 6 7 , cos 9 1 , cos 7 ) = (C+ - cos 7 ) (cos # 7 - C~) , (11) 
= cos (9 7 =f # 7 ) = cos6 7 cos# 7 ± (l — cos 2 Q-^ 1 ' : - 11 
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and 7 is the angle of Compton scattering from p 7 to p 7 , 

cosG 7 (p 7 ,p 7 ) = 1- ^ + ^ (12) 

/>7 i^ 7 

(here p 7 < p'). In (fTUj) . we keep 7 , # 7 for easier integration over the angular variables, but before 
appplying d/dr in the LHS of ©, cos# 7 should be expressed in terms of p using ([7]). 

The last term in the RHS of ([9]) describes the total kinetic departure of p 7 -photons due to 
their scattering with the total collision frequency, v(p-y = m e ca), proportional to the total K-N 
cross-section [4], 



2cn m Z air c7KN 

2 + 8a + 9a 2 + a 3 (2 + 2a - a 2 ) In (1 + 2a) 



47rn m rocZ air 



a 2 (1 + 2a) 2 2a 3 



(13) 



Equation Q should be supplemented by the boundary conditions. In accord with (|8b|). at a 
small distance r that satisfies a <C r -C Z = c/z/, the photon distribution function radiated by the 
point-like radioactive source can be written as 

k 

A + (^7> P) | r _> = 5^ (PV P) 1^0 ' 
s=l 

f + ( xi _ 4 0) ^ ( cos °i - l ) 6 (Pi ~ Ps ) - {p} 6 &i ~ Ps ) fi^n 

/ 7S ^7»P;| r ^ 87T 2 Cr 2 p2 87T 2 Cpp2 ' W 

where 1^ = 3.7 x 10 10 s _1 Q (Ci) is the total number of radioactive decays per unit time, X)s=l 
denotes summation over all kinds of primary photons radiated during one disintegration of the 
radioactive material (usually, k = 1,2,3) and 6 (p 7 — p s ) expresses the fact that all radiated photons 
have discrete energies £ s = cp s (e.g., for 27C0: k = 2, S% = 1.332 MeV and £2 = 1.173 MeV). 
Equation (fT4"|) represents the inner boundary condition at r — > applicable to the outgoing photons. 
For the ingoing photons, we use the natural condition of /~ (p 7 ,p) — > at r — > 00. Between 
consecutive collisions, some ingoing photons pass the minimum distance r = p and become outgoing. 
This yields the matching condition between the two groups: 

fj (r,P-y,p = r) = f~ (r,p 1 ,p = r). (15) 
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Further, all collisional processes lead only to decreasing photon energies, so there are no photons 
with energies exceeding the initial ones, /+ (p > max(p s )) = 0. Lastly, we assume that at p = 
(cos# 7 = ±1) all momentum distributions behave regularly. 

The RHS of Q satisfies automtically the condition p 2 dp-y ^(Lj-f^—vf!^ 1 ) (r 2 — p 2 ) pdp = 
0, which provides in accord with (I8b[) conservation of the total photon flux through any surface 
surrounding the radiation source. This flux conservation is due to the fact that kinetic equation 
([9]) includes only photon scattering. Photoelectric absorption and other possible photon losses can 
also be included in ^(j> 7 ) by a simple addition to the RHS of (|13j) . Those effects would break the 
total flux conservation. 

Kinetic equation is solved analytically in Appendix [Bj The complete solution is given by an 
infinite series of multiply scattered photons (|B.2p . where /q" = (/^~) prim is given by (jB.ip . f^ >1 is 

determined by recurrence relations (|B.4p with the linear integral operator L r defined by (jlOp . This 
general solution, however, involves an increasing number of nested integrations which for n > 2 
are difficult to calculate. Below we restrict our treatment to distances r from the radiation source 
located within the mean free path for 1.2 MeV 7-quanta, r < l s = c/v(p s ) ~ 120 m (the "near 
zone"). At these distances, one might expect two dominant photon groups: the primary 7-quanta 
beams (n = 0) and singly scattered photons (n = 1). The momentum distribution for singly 
scattered photons is obtained in Appendix IB.2I However, a mere cutoff of the infinite series (]B.2p 
to its first two terms, f$ and f x given by (IB.17j) - (IB.22p . becomes less satisfactory at distances 
r ~ l s because further scattering of singly scattered photons breaks the total flux conservation. 
Below we modify the two-group description of photons to make it exactly flux conserving. 



4 Photon energy distribution in the near zone 

To calculate the total electron production rate dN e /dt given by ([6]), we need to determine the 
spatially dependent photon energy distribution at various distances from the radioactive source, 
-F 7 (a), where a = p^/{m e c) is the normalized energy. In this section, however, it is more convenient 
to use 

J ' ± Jo (r z — p A ) ' 

which differs from by the normalization factor following from dn^ = F(p~)dp™ = F{a)da: F 1 {a) = 
m e cF(p^). In terms of F(p^), ([6]) becomes 

f r\N \ 1 r°° 

I ) ~ 2.09 x 10 10 s- 1 / Pl F( Pl )d Pl . (17) 

V dt y prim m e cJo 

4.1 Primary photons 

For the primary photons, we have the solution (jB.ip valid for arbitrary distances from the radiation 
source. The exponential factors e~ r ^ ls there describe the beam attenuation due to Compton scat- 
tering. For the corresponding F{p 1 ), photon density (^7) Drim , total flux (-£y) Dr i m and the electron 
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production rate (dN e /dt) im , we obtain simple formulas 

s=l s=l 



^) » 2.09 x 10 10 s^ 1 V -EL. e -r/ls 

dt J prim 4vrcH ^ m e c 



w 4.017 x 10 5 s" 1 cm" 3 V^(MeV)e- r/is . (18) 

[ r ( m )l s= i 

4.2 Scattered photons: extended n — 1 solution 

Each primary beam of radiated photons denoted by s produces its own set of scattered photons 
characterized by the corresponding partial distribution function /+. Bearing in mind that the 

distribution of all photons is given by = Y^, s =i fys an d 



7s 

F iPi) = 2^ F s(Py), Fs(p 7 ) = —— 2^ / — T 1777^- 

«=l r ± - 70 ( ri - P 2 ) 



(19) 



below we consider separately the partial distributions fj~ 3 and the corresponding partial moments. 

As a primary beam of photons leaves the radiation source it starts producing scattered photons 
with the rate, which in the close proximity to the radiation source, r <ti l s , grows linearly with dis- 
tance. Convolved with the 1/r 2 spherical divergence, the population of scattered photons becomes 
proportional to 1/r. Our simulations show that this spatial dependence for the energy distribution 
of scattered photons roughly holds through the entire near zone, r < l s . Until the low-energy 
photoelectric absorption becomes important, a partial photon flux of all related photons remains 
constant, J~ 8 = Iy . Each partial primary-photon beam s decays exponentially due to Compton 
scattering, (^7s) prim = I^e~ r / ls , so that at r > 0.7l s the corresponding flux of scattered photons, 

Zy ^ (l — e~ r / ls ), dominates. The corresponding photon density, (n 7 ) = I7 (l — e~ r l ls ) /(v rs ), 
where (v rs ) is the mean radial velocity of scattered photons, starts dominating noticeably closer to 
the radiation source due to the more isotropic moment um distribution ((tVs)scat 

< c) compared to 

that of the almost perfectly directional primary beam ((vrs)prim ~ c). The same pertains to the 
electron production rate, albeit the distance where dN e /dt produced by scattered photons starts 
dominating at a larger distance r than that for the photon density due to the weighting factor p 7 in 
the integrand of (|17p . (This factor is biased to higher energies where scattered photons momentum 
distribution is more directional.) Thus electron production due to scattered photons may become 
important deeply in the near zone, r <l s . 

We simplify the treatment of scattered photons in the near zone as follows. We will consider all 
scattered photons, n > 1, as one unified group and extend for this group the singly-scattered (n = 1) 
photon momentum distribution given by (|B.17|) and (|B.21|) with one simple modification. In the 
exponents of the attenuation factors (1B.18|) and (|B.22h . we will retain the terms r s /l s responsible 
for the primary beam attenuation, but will ignore [±(r 2 — p 2 ) 1 ^ 2 ± (r 2 — p 2 ) 1 ^ 2 ]/! (p 7 )- The reason 
is that the latter terms describe the loss of n = 1 photons due to their further scattering. For the 
entire unified group of all scattered photons, this represents just a particle redistribution and hence 
must be ignored. This will provide automatically the total photon flux conservation. 
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The extension of the n = 1 solution to all scattered photons will introduce two inaccuracies. On 
the one hand, multiply scattered photons decrease their energy as n increases, so that photons with 
n > 2 leak into the lower-energy range £ < £^J n = £ s [l + 2p s /(m e c)] ~ 0.2 MeV, as discussed in 
the end of Appendix section IB. 11 This reduces the photon population in the regular n = 1 energy 
range, £^ n < £ < £ s - Due to the weighting factor p 7 in the integrand of (fT7|) . this inaccuracy 
in the extended n = 1 solution would lead to overestimated dN e /dt. On the other hand, multiply 
scattered photons have a more isotropic angular distribution in the momentum space that leads 
to the reduced average photon radial velocity and, hence, to the increased photon density. As a 
result, the more directional extended n = 1 solution would, on the contrary, underestimate dN e /dt. 
We expect the two inaccuracies to partially balance each other. As discussed below by comparing 
our analytical results with computer simulations, our simplified analytical model describes well the 
spatial dependence of the photon energy distribution within £^ n < £ < £ a , practically in the entire 
near zone (r < 100 m). Thus it can be employed for reliable estimates of the electron production 
rate in the near zone. 

The extended n = 1 solution is even simpler than the original one because (IB. 18ft and (IB. 22ft 
reduce to a common attenuation factor, A ext = ^4^ t = exp (—r s /l s ). One can easily check that the 
modified solution conserves exactly the total photon flux through any closed surface surrounding 
the radiation source. Indeed, the contributions of the ingoing and outgoing photons at r < r s (see 
Fig. H]) cancel each other. The remaining partial fluxes of the outgoing photons in the outer sphere 
yield 1^ (l — e~ r / ls ) , in full agreement with the flux conservation. Unfortunately, the contributions 
from the ingoing and outgoing photons within r < r s to the photon energy distribution and dN e /dt 
have no opposite signs and add rather than cancel. This leads to more complicated calculations 
outlined below. 

Using the modified n = 1 solution, we calculate now the energy distribution of scattered photons. 
According to the above, the partial momentum distribution assoiciated with a beam s is 



v r(°) 2 / ( \ 

Zj a \ v ij n m r^m e r s r ^ j r s [p,p*y,p s ) 

Airp 2 p 2 p 



Jjs- A„„2„2„2 exp ( - , ], (ZUJ 



where 



ife) = ^)- (21> 



and r s (p,p^,p s ) = p/L 1 / 2 , introduced in (IB.14j) . describes the distance where the scattering of 
primary photons took place. Introducing L(p^,p s ) = sin 2 7 , as in (|B.9p . we obtain 

fi(, 7 ) = ^y° v (?i + el- l )! 5° r * ^^ = it^) (22) 

2rp2 \ Ps Pj J { S hack if Pl < p cvs = 1+p ff (meC) 

where the boundary p 7 = p CTS corresponds to L = 1 (6 7 = 90°) and 

1 f rLl/2 exp(-r s /l s )dp 



VLJo {r 2 -p 2 ) 1/2 ' 




11/2 exp(-r s /l s )dp + 2 r exp(-r s /l s )dp \ ^ 



(r 2 -p 2 ) 1/2 Jrtf-f* (r 2 -p 2 ) 1/2 

The function Sf OT describes photons scattered forward (0 7 < 90°). It includes only outgoing photons 
with r > r s (0 < p < rL 1 / 2 ). The function 5b ac k describes photons scattered backward (0 7 > 90°). 
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In the "inner" zone, r < r s (rL 1 / 2 < p < r), it includes both ingoing and outgoing photons 
with equal contributions, while in the "outer" zone, r > r s , it includes only outgoing photons. 
The integrals over p describe the original scattering of photons from the primary beams, where 
the distance of scattering, r s = p/L 1 / 2 , depends on the angular parameter p and on the photon 
momenta via L(p 7 ,p s ). In the near zone, r/l s < 1, all outer photons at given r were scattered 
at distances r s < r. The situation for the inner photons is different. Some ingoing photons at 
given r could be originally backscattered at large distances, even r s 3> r (for L < 1), well beyond 
the near zone. The total contribution of these photons, however, should be relatively small due 
to exponential attenuation of the primary beams described by exp(— r s /l s ) = exp[— p/ (L l / 2 l s )}. 
Temporarily introducing 



P(L,Xs 



Li/2 



exp (— X s t) dt 



2\V 2 



X s (r,p~ p p s ) 



IM 2 L 



m e c m e c 



Pi 



Ps 



-1/2 



Ps 



m e c 
Pj 



+ 2 



-1/2 



and t = p/r, we can rewrite the functions Sf or and S'back as 



Sfo 



P(L,Xs) 



S'back 



2P(1,X S )-P(L,X S 



(24) 



(25) 



(26) 



The transition from Sf or to S'back at L = 1 corresponding to p 7 = p crs is smooth. For example, 
P(L,0) = arcsin(L 1 / 2 ), while 2P(1,0) - P(L,0) = vr - arcsin^ 1 / 2 ) is the continuation of the 
multivalued function arcsin(L 1 / 2 ) from its first quadrant (0 < arcsin(L 1 / 2 ) < vr/2) to the second 
one (n/2 < arcsin(L 1 / 2 ) < it). A similar smooth transition takes place for non-zero A s . 

If A s < 1 then the scatter distances r s lie within the near zone. However, as mentioned above, 
ingoing photons scattered backward can arrive from arbitrarily large distances r s , so we should 
consider any possible values of A s > 0. All such photons are described in (|26j) by the function 
P (1, X s ) shown in Fig. [5] and expressed analytically as 



7T 



P(1,A S ) = -(J (A S )-L Q (A S )) 
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(27) 



where Iq(x) and L${x) are the modified Bessel and Struve functions [10J. For a simpler analysis, 
one can approximate P (1, A s ) in terms of elementary functions by combining the power series (|27l) 
to the fifth-order term with the large-A s asymptotics, P (1, A s ) — > l/X s , as 



P(1,X S 



1 + 



-A, 
1/A 



1 I A s I A s 

1 i" 9 ^225 



if 
if 



A s < 2 
A, > 2 



(28) 



This piece-wise approximation mimics the actual P(1,A S ) almost ideally except for a minor kink 
at 1.5 < A s < 4, as seen in Fig. [5j 

Unlike P (1, A s ), the function P(L,X S ) in (|26l) describes only outgoing photons scattered at 
distances r s closer to the source than the given distance r. According to the near-zone restriction, 
r/l s < 1, in the entire integration domain of (j24[) . X s t < 1 holds. This allows us to expand the 
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Figure 5: Function P(l, A s ) given by (f27l) (solid black), along with its approximation given by ([28]) 
(dashed red). 
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Figure 6: Photon energy distribution for unshielded 27 Co at various distances r. The analytical 
solutions are shown by solid lines. 



exponential function in the integrand of P (L, A s ) in powers of (\ s t) to the fifth power (applicable 
for \ s t < 2). This yields 



P(L,A S )~ (l + y) arcsinv / L-A s (l + ^ + ^ (l - VT=L 

A 2 (4 + 3L) 3\ S \/~L 
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(29) 



Returning from the temporary notation A s = r /{l s ^/~L) (125p back to the distance r, we obtain 
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(30) 



(31) 



where, according to (I27p - (128p . the first term in the RHS of (131 f) can be approximated by elementary 
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MCNPX Simulation CO-60 in air (1 cm thick Pb shielding @ 99 cm) 
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Figure 7: Photon energy distribution corresponding to Fig. [6l but for the 1 cm thick spherical lead 
shielding placed at 99 cm from the source. 
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functions as 
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i s VlJ \i s Vl 



TT 

7z 

J^( 1 +|l) 2 -S( 1 +9& + ^) if VL>r/(2l s ) 

I 2l s /r if VL<r/(2l s ) 

The function Sf OT describes photons scattered forward (0° < 7 < 90°, p CTS < p 7 < p 

max — Psji 

while Sb a ck describes photons scattered backward (90° < 7 < 180°, p mm = p s /[i- + 2p s /(m e c)] < 
P'y < Pcrs)- The two functions are smoothly connected at L = 1 (0 7 = 90°, p 7 = p C rs)- In spite of 
singular denominators L~ q in ([30 ]) - ([32]) . 5f or and S^ack behave regularly at L = (0 7 = 0°, 180°), 
except for r = 0. At r <C 5f or (0 7 = 0°, p 7 = p s ) approaches 1, while S^ack (©7 = 180°, 
P~f = Pmm) increases with decreasing r as 2l s /r. This 1/r divergence of Sb ac k(-£ = 0) at r — > is 
due to accumulation of photons scattered almost strictly backward from distances r s ^> r where the 
primary beams still remain largely unattenuated, r s <C l s - As r increases, even remaining within 
the mean free path l s from the radiation source, the exponential attenuation of the primary beams 
at larger distances r s > l s makes such accumulation less significant. 

Inserting ([30]) and (j3T|) into (j22|) . with L given in terms of p 7 ,p s by (|B.9|) . i.e., 

/ \ 1/2 / \ 1/2 

/— 7 / m e c m e c \ ' I m e c m e c x ' 
VL = sin 9 7 = —2 —2 + 2 

V P~f Ps J \ Ps P-y 



vl — L = Icos © n 



1 — H — 



(33) 



we obtain the complete approximate analytical expression for F{p~) (|19p in the near zone. 

Figure [6] shows photon energy distributions for unshielded 27 Co calculated by the MCNPX 
code (URL: http://mcnpx.lanl.gov/) and the analytical solution given by ([22]) and ([30]) - ([33]) for 
several distances r (shown near the curves). The two long spikes on the right correspond to 
the two monochromatic primary photon beams with S\ = 1.332 and £2 = 1.173 MeV, described 
analytically by (jB.ip . We see an excellent agreement between the simulations and theory in the 
entire energy range of singly scattered photons, £2/ (1 + 2£2/m e ) « 0.21 MeV < £ < 1.332 MeV. 
At £ < 0.21 MeV, the simulations show lower-energy multiply scattered photons excluded from the 
extended n = 1 model. In the near zone, especially due to the weighting factor p 7 in the integrand 
of (|17p . the lower-energy range makes a small relative contribution to dN e /dt. 

Figure [7] shows simulated photon energy distributions for the same radioactive source but 
shielded with a thin spherical lead cover. The clearly pronounced broad bumps at the high-energy 
tails of the photon distribution, especially at the closest proximity to the source (r = 1 m) are due 
to primary photons scattered within the shield material and leaked outside. This scattering could 
also be described by the general solution (|B.ip - (IB.4j) with L r defined by fllOh . where the collisional 
mean free path for photons in the shielding metal is many orders of magnitude smaller than that in 
air. Unfortunately, for dense shielding the near-zone approximation is hardly applicable. Besides, 
photoemission absorption may play a role. This makes analytical calculation much more difficult. 
Comparing Figs. [6] and El we see that the effect of shielding on the energy distribution is especially 
pronounced in a close proximity to the source. For distances > 50 m, the photon energy distribution 
in both shielded and unshielded cases look similar. 
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4.3 Electron production rate 



Now we calculate the electron production rate (fT7|) due to scattered photons with the energy 
distribution function given by (|22p . or 
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where a s = p s / (m e c), L = 1 — (1 — l/a+ l/a s ) 2 and 5f or , S'back are given by ([30|) -([32 |) . We restrict 
calculations to the radiation material 27C0 with E\ = 1.332 MeV and E2 = 1.173 MeV. 

When calclulating the production rate for scattered photons, (dN e /dt) sca _ t , it is convenient to 
compare it with production rate for primary photons given by (|18p . Straightforward numerical 
integration shows that for each s = 1,2 a simple approximation 
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(34) 



holds to a better than 5% accuracy. To compare with {dN e / dt) pr \ m , we express l s = c/v s , with v s 
in terms of a s by (fl~3j) , and obtain 
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As a result, the total electron production rate becomes 



dNe 
dt 
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(35) 



where the last term in the RHS brackets represents the ratio between the production rate due to 
the scattered photons to that due to the primary photons (for each primary beam s) and 
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(36) 



Equation (|35p shows that the scattered-photon contribition to the total free-electron production 
rate decays with distance as e -°- 42r /^ i.e. noticeably slower than that of the primary photons, 
e -r/h/ r 2, For 60 CO) we have i x ra 127 mj as „ 2.607, Ml ~ 0.72, and Z 2 ^ 119 m, a 2 ~ 2.295, 
Hi ~ 0.77. The corresponding (dN e /dt) I£kd is shown in Fig. [HJ along with the primary photon 
rate given by (|18p . This Figure shows that the latter is a good approximation for r < 40m. At 
larger distances, dN e /dt due to scattered photons becomes noticeable and dominates starting from 
r ~ 100 m. For simple estimates, we can extrapolate equation (|35p to larger distances, but for our 
application only moderate distances, r < 150 m, matter (see below). 

Figure [8] also shows that a moderate amount (~ 1 mg) of a concealed radioactive material like 
27C0 can provide the total rate of electron production well above the background level (Q ra d ~ 20 
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pairs s _1 cnr 3 ). The enhancement factor a ra d = {dN e /dt) raA /Q ra ^ exceeds unity within 150 m 
from the radiation source. This distance is crucial for implementation of our technique because 
... As discussed in Ref. [3], all free electrons will eventually end up forming negative ions, but the 
collisional detachment of some electrons for a brief time make it sufficient to initiate air breakdown 
within the focal spot of the THz gyrotron during pulse. Equation ([35]) provides a reliable estimate 
of the electron production rate at least within the distance of 150 m critical for implementation of 
our method. 

Conclusion: In this paper, we analyze the production of free electrons by 7-rays leaking from 
radioactive materials. This is important for implementation of a recently proposed concept to 
remotely detect concealed radioactive materials pQ, [3]. We have performed our analysis specifically 
for the most probable radioactive material 27C0, but the general approach can be applied to any 
7-ray sources. In our theory, we have included free-electron production by both primary 7-quanta 
radiated by the source and penetrated through the container walls and those scattered in air. 
Equation (f35|) and Fig. [8] give a reliable estimate of the total production rate (dN e /dt) rSbd /Q Ta ,d and 
the corresponding radiation enhancement factor, at least within 150 m from the radiation source. 
At these distances, the scattered photon contribution to (dN e /dt) vad is comparable or even exceeds 
that of the primary photons. Even for amount of radiation material as small as ~ 1 mg, the 
total production rate may exceed significantly the ambient ionization rate. During the gyrotron 
pulse of about 10 microsecond length, such electrons may seed the electric breakdown and create 
sufficiently dense plasma at the focal region to be detected as an unambiguous effect of the concealed 
radioactive material. 

Acknowledgement: This work is supported by the US Office of Naval Research. 



Appendix 

A Deriving kinetic equation ([9]) 

The Boltzmann kinetic equation for freely propagating unpolarized photons colliding with air 
molecules contains two parts: the collisionless part describing the undisturbed propagation of 
photons between two consecutive acts of scattering or absorption and the collisional operator [8]. 
In the assumed spherically symmetric approximation, the stationary kinetic equation for fj^ (r,p, p) 
with p given by ([7]) can be written as 

± J 1 p 2 V 12 df±{r lP ,p) 
\ r 2 J dr 

The LHS of (|A.ip contains no p, /9-derivatives because these variables are invariants of the collision- 
less motion. The collisional term [d/^F/di] describes scattering, absorption and other possible 
effects responsible for collisional changes in the photon momentum distribution. In the general 
case, it has a complicated integral form which can be simplified in the spherically symmetric case. 

In this appendix, we obtain the explicit simplified form for [df^/di\ We will start directly 
from the original probabilistic expressions like [8j. Since the photon- molecule collision is a strictly 



dt 



mil 



(A.l) 
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Distance from radiation source (in meters) 



Figure 8: Free-electron production rate for 1 Ci (0.9 mg) of 27C0 as a function of the distance from 
the raioactive source. The red dashed line shows the primary-photon rate given by (|18p . while he 
black sold line shows the total rate given by (|35j) . 
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local process where the photon transport plays no role, in the RHS of (|A.1|) we can temporarily 
return from p back to the momentum polar angle # 7 . 



tf 7 )d 3 p' 7 d 3 p e , 

(A.2) 

p 7 )c? 3 ^(iVe- 

Here Idep describes the total collisional departure of photons from the kinetic volume d 3 p y around 
given p 7 with scattered photon momenta in the domain d 3 p'^ and releasing electrons in d?p' e around 
f? and p e , respectively. In the general case, the departure term, hep-, can include both photon 
scattering and absorption. The term I arr describes the arrival of photons and released electrons into 
the volumes d 3 p^ and d?p e around p 7 and p e , respectively, if the incident photons before scattering 
were at d 3 p 7 around p^. Since we neglect any collisional photon production like Bremsstrahlung, 
fluorescence, etc., I arr includes only photon scattering. With no photon losses included in Idep, 
(|A.2|) provides automatically the particle conservation, J (df-y/dt) co nd 3 py = 0. 

The function W of two arguments p f and p 7 represents the probability per unit time of the 
collisional process. This function is proportional to the air molecular density, n m = n a /2. In the 
general case, the two arguments of W are not interchangeable. To identify the correct order of the 
arguments, the time sequence of the collision is denoted by the arrow. The two delta-functions in 
the integrands describe the energy and momentum conservation. Since we assume energies well 
above the K-shell edge, the arguments of the delta functions include no small energy spent on the 
molecule excitation and ionization. They also neglect the kinetic energy and momentum of air 
molecules. As explained in section O we presume Compton scattering in the K-N approximation. 

To obtain the explicit expression for W, we relate it to the corresponding differential collision 
cross section, da/dQ'^, given for the K-N theory by Eq. ([2]). To find the expression is easier through 
the departure term Idep because, unlike I arr , it does not involve the photon distribution function 
in the integrand. Integrating (|A.2j) . first, over f? e with eliminating 5(p r ^ f + p' e — p-y) and then over p'^ 
in d 3 p^ = p'^dp'^dtl'^ with eliminating the remaining delta-function, we obtain 

2cZ ^ m ^=W { 0^ *,) \ £ , Jc ^ +m ^ (l + ~J§j\r (A.3) 

Here 2Z a ; ir n m is the total density of bound electrons and d£' e jdp'^ should be calculated from the 
momentum conservation, p 7 = + p e , but before applying the energy conservation, p 7 + m e c = 
£' e /c + p' T As a result, we arrive at 

„, / — / \ dep Z^rimr^meC 3 (p-y Pj 2 \ 

w ^^^)= e*p^ U + ^" 1 + cos Q V' {AA) 

where 7 is the scatter angle from p 7 to given by ([3]). This results in Idep = ^(p^fyfr,^), 
where v{py) is given by (fT3l) . 
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Now we proceed to calculating the more complicated arrival term, I arr . Switching in (|A.4|) the 
arguments of W, we obtain 
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arr ^air n m r \' m e ( ? 
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(A.5) 



where cos 7 in terms of p-y and p' > p 7 is given by (fl2j) . After inserting (|A.5|) to (|A.2[) . as above, 



the first step to calculating J arr is to integrate over p e with eliminating 5 (p 7 + p e — p 7 
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This yields 
(A.6) 
within the 



Unlike Idcpi however, the arrival term I arr contains the distribution function /' v . . .. , 
integrand. This does not allow us to easily integrate over dp'^ in d 3 p'^ = pfidpLd£l' in order to 
eliminate the remaining delta function. However, the fact that the distribution function f depends 
on p 7 and # 7 , rather than on the axial angle 7 suggests eliminating the remaining delta function 
via the integration over (ft' using a simple geometric relation 
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Given p 7 , # 7 ,p 7 , $ 7 , this yields E 12 (</> 7 ). Integration of in (jA.2j) over (ft 1 using 
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The factor 2 in ()A.9p originates from existence of two symmetric values of (ft'^ 
argument of the delta-function to zero (the energy conservation) 

E e = p' + m e -p 7 . 



(A.7) 

(A.8) 

(A.9) 
that set the 
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The specific value of <ft' c in (|A.8|) is inconsequential. The differentiation in the normalization fac- 
tor \dE e / dcft'^ should be done using ()A.6|) and (|A.7|) . but before applying (jA.lOp . This yields 
E' e dE' e /d(ft' = p'^Pj sin 9'^ sin # 7 sin (eft' — 7 ). Using (|A.7|) again in order to express sin {(ft' — 7 ) 
and then, applying the energy conservation relation (|A.10|) . we obtain \dE' e / dfi^ 1 = E' e j fap^S 1 / 2 ), 
where 



S (cos 7 , cos # 7 , cos 9j) 

= 1 — cos 2 6-y — cos 2 0' — cos 2 9j + 2 cos 7 cos 9' cos 6> 7 (A. 11) 

Implementing the entire procedure, we arrive at the integro-differential equation ([9]), where S is 
recast as (jlip . This fully symmetric function of three cosine arguments has a simple physical 
meaning. The restriction S > covers all possible directions of the incident and scattered photons 
characterized by the pairs of polar and axial angles 7 ,<3? 7 , # 7 ,</> 7 , # 7 ,</> 7 , which are coupled by 
geometric relations like (|A.7p . The photon energy losses and the scattering angle 7 are rigidly 
coupled by (fT2"jh A photon with given 9' p' can be scattered into another photon with given 6* 7 , 
Py, provided the photon energies are within the allowed domains. This freedom is reached by a 
unique choice of two symmetric values of (</> 7 — (ft'). 
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B Solving kinetic equation ([9]) 



B.l General solution of multi-scattered photons 

The LHS of kinetic equation ([9]) is simple, but its RHS is rather complicated, as described by (fTOj) . 
Our approach to solving the kinetic equation is based on the physics of the process. The primary 
photons with the monochromatic momentum distribution approximated by (|14p start efficiently 
scattering by air molecules within a distance of r < / = cjv from the source. For photon energies 
between 0.3 and 1.3 MeV, the maximum K-N differential cross-section is reached near the forward 
scattering, 7 = and gradually decreases to the backscatter, 7 = n, as shown in Fig. [2j This 
means that photons scatter mainly through large angles, rather than broaden their delta-function 
angular distribution around the primary beam. This suggests considering the primary photon beam 
as a separate group that continues keeping its delta-function distribution over long distances from 
the radiation source. For the primary photons, we can neglect the first ("arrival") term in the RHS 
of (J9|) but must keep the last ("departure") term. Using the boundary condition (Till) , we readily 
obtain for the primary photons the exponentially decaying solution, 
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This solution gives the entire local source for scattered photons. The solution for all propagating 
photons, including the primary beam, can be sought as an infinite series 
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where /q~ = (/^") rim and the subscript n denotes photons scattered from the primary beam n 
times. For each individual group of n-scattered photons, we obtain a recursive kinetic equation 
where the source for scattering is provided by the (n — l)-th group, 
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where L r is the linear integral operator defined by (|lUp . The infinite series given by (|B.2p converges. 
If v includes no absorption then Eq. (|B.3P provides automatically the constancy of the total photon 
flux through any closed surface surrounding the radiation source. Using the boundary and matching 
conditions for the outgoing (+) and ingoing (— ) particles described above, we obtain the exact 
recursive solution for in terms of f^_ 1 , I = c/u, 
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This solution holds for p < r, while for p > r all f£ = 0. As mentioned above, the integrations 
for L y in (jlpp can be performed it terms of (cos 9' J , but before inserting the results to (IB. 4ft they 
have to be transformed into functions of p using (|7|). In the integrands of (|B.4|) . the distance r is 
replaced by an internal radial variable y in order to avoid a confusion with r outside the integral. 
The exponential factors describe (n — l)-photon losses due to their scattering into the n-group with 
|(r 2 — p 2 ) 1 / 2 ± (y 2 — p 2 ) 1 ^ 2 \ describing the distance between two points on the r- and y-shells along 
a straight-line photon trajectory corresponding to given p (the signs ± depend on whether the two 
points are on the same or on different sides with respect to the trajectory median r = p, as shown 
in Fig. |3J). The summations over the ingoing and outgoing particles in ()B.4p take into account the 
fact that photons can be scattered through arbitrary scatter angles, so that any ingoing or outgoing 
particle can be scattered into each of these groups. 

The explicit solution for scattered photons with arbitrarily large n can be found by consecutively 
applying recursive relation (IB.4j) starting from n = 1 with = (/^) prim given by Eq. (IB.lj) . At 
any given distance r, photons with sufficiently large n are produced mainly by ingoing photons that 
originate from the primary source beam at much larger distances and only accidentally come back. 
One might expect that they will come back less frequently as n increases, so that the infinite series 
(|B.2p will converge. Additionally, as n increases, the photon energy distribution gradually moves 
to lower energies. This effective "cooling" of scattered photons will eventually put high-n scattered 
photons into the keV range where the photoemission absorption of photons starts dominating. 

We can easily estimate the energy spread for n times scattered photons. The maximum energy 
of scattered photons, p 7 = p' — I, where I is the ionization potential, I ~ 15 eV <gC p s , is reached 
where # 7 is close to 0°, but not too close [U Sect. 2.6.2.1]. As a result, after n scatterings of 
a primary 7-quanta with p 7 = p s , the maximum scattered photon energy reduces only weakly, 

(n) (n) 

Pmlx = Ps — n I) while in the framework of the K-N theory pmax = Ps- On the other hand, the 
minimum energy of scattered photons is reached for strict backscatter, # 7 = 180° jlj. After n 
scatterings, the minimum energy, p m ^ n , decreases by the following rule, 

(n-l) 

starting from p^ m = Po- Temporarily denoting Y n = m e / '(2p^ n ), we obtain from ()B.5|) Y n = 
1 + y n -i, so that Y n = n + Yq, or 

pW = * = ^ . (B.6) 

lmm l + 2np s /(m e c) 2n(l + m e c/(2n Ps )) v ; 

The energy of n-times scattered photons spreads between p^ n and pmix- This spread gradually 
increases with n due to lowering p^ n - For p s ~ 1 MeV, 2p s / (m e c) ~ 4, so that after 2-3 scatterings 

the minimum photon energy p~- n becomes almost inversely proportional to n, p^ n — m e c/(2n) ~ 
0.255 MeV jn. Hence, after ~ 10 scatterings, the minimum energy falls well into the range where 
photoemission absorption starts dominating. This low-energy photon absorption creates an effective 
sink for multiply scattered photons. 



B.2 Solution for n = 1 

Here we consider the lowest-order scattered photons, n = 1. In the relatively near zone, r < l(p^), 
the cutoff of the infinite chain by the first term of the series represents a good approximation for 
the entire photon distribution. 
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For n = 1, we obtain from (|B.4p . 
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Here (4rr)i is given by 
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where cos 7 (p 7 ,p 7 ) is given by (fT2|) and, accordingly, 
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In accord with (IB.ip . we obtain 
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Unlike general equation (IB.4j) . (IB.7D contains no summations over ± because the primary photon 
beams include no ingoing photons. Integration in (|B.8|) over cos # 7 should be done with care 
because the straightforward elimination of <5(cos6> 7 — 1) with S (cos 9 7 , cos 6> 7 , cos 9' 7 = l) leads to 
a singularity. To avoid it, we apply the following. Regardless of the specific values of cos 7 and 

cos6> 7 , the integral J* c _ (l/\^S)d (cos# 7 ) equals n, so that we obtain 
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where l s = l(p s ) = c/v(p s ), 

, \ P f m e c m e c\~ 1/2 ( m e c m e c V 1/2 

r s (p,p 7 , Ps ) = — = p^—-—) (___ +2 j . (B.14) 

These expressions have a simple physical meaning. Photons from the radiated monochromatic and 
infinitely narrow beams are scattered into angles rigidly coupled to the energy of the scattered 
photon. Given both the energy and p of the scattered photon, one can find the distance where such 
act of scattering took place, r = r s . The delta- function 5 (r — r s ) demonstrates the uniqueness of 
this scattering distance. The factor e~ Tc l la describes the total attenuation of the primary beam s 
on its way from the radiation source to r s . 

Now we insert the expression for (I a rr)i into (|B.7|) . integrate with eliminating 5 {y — r s ) and use 
relation (|B.14[I : y r s , {y 2 — p 2 ) 1 ^ 2 —> (r 2 — p 2 ) 1 ^ 2 . The result depends upon whether the photon 
momentum p 7 is less or greater than a critical value, p crs , corresponding to the 90°-scattering 
(cos G 7 = 0) from the primary beam s. For 

Pcrs = — % r < P 7 < Ps (B.15) 

1 +p s /(m e c) 

all singly-scattered particles are outgoing particles caused by forward scattering (0 7 < 90°), while 
for p 7 < p CTS scattered particles are either ingoing or outgoing particles caused by backscatter 
(90° < 6 7 < 180°). 

Depending on whether r > r s or r < r s , we will separate between 'outer' and 'inner' photons. 
For p 7 > p crs , all scattered photons are outer (r > r s ) outgoing photons, 

/ \ 1/2 / \ 1/2 

T i/9 / f^eC m e c\ ' I m e c m e c \ ' .„ 

p<L 1/2 r= — — — — + 2 r, B.16) 

V Pi Ps J \ Ps P-y ) 

originating from only one scattering location with given r s . We have then 

= ^I^n r 2 m e j^(p_ x + P_ ± _ \ 

8-rrp 2 p 2 \Ps p-y J pi 

where the attenuation factor is 

A+. (r,p 7 ,p) = exp {[(r 2 - p 2 f 2 - (r 2 - p 2 f /2 )/l (p 7 ) - r s /l s ) (B.18) 
For p 7 < p CTS , or more precisely, for 

Pmin = -77 r < P 7 < Pcrs = ~ % r (B.19) 

l + 2p s /(m e c) l+p s /{m e c) 

there are both outer outgoing and inner (r < r s ) photons, 

L x l 2 r < p<r, (B.20) 

either outgoing or ingoing from two different locations of given r s . For the outer photons defined 
by (|B.16p . we have the same solution as (|B. 17|) . For the inner photons, we have two solutions for 
the either ingoing or outgoing photons, 

97 t(0) 2 fc / \ A± 

e± __ ^^air-Lj n m r Q m e ST~^ I Pi . Ps T , \ \ ^int^s 



87rp2p2 y ps p ^ J Pi 
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which differ by the attenuation factors only, 



Af nt (r,Py, P) = exp [- (r 2 - p 2 f 2 /I (p 7 ) T (r 2 - p 2 f 2 /I (p 7 ) - r s /Z s ] (B.22) 

1 /2 

There are different signs in front of v (r 2 — p 2 ) /c due to the fact that for the outgoing 
photons the two points corresponding to r and r s lie on the two opposite sides from r = p, so 
that the two distances (r 2 — p 2 ) 1 ^ 2 and (r 2 — p 2 ) 1 ^ 2 add, rather than subtract, as they do for the 
ingoing photons. 

Note that the only r-dependence in the above distributions is due to the attenuation factors. As 
functions of the angular variable p = L 1//2 r, contains only a singular factor r s /p 2 oc 1/p. Since 
all integrals for calculating the photon densities and fluxes contain d (cos 7 ) oc pdp this singularity 
is integrable. 

For n > 1, the formal general solution (|B.2j) with the recurrence relations (|B.4j) involves an 
increasing number of nested integrations and becomes complicated. The solution for n = 1, given 
by equations fBTf|) - ([Bl22]> . will be the basis of our analysis of the photon propagation in the 
near-source zone. 
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